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in the case of ordinary bosons, p = 1. Two non- negative real parameters arise because "a6" 
and "5a" are fundamentally distinct pair operators of charge A-quanta and charge '—1', 
^1 B-quanta parabosons. hi 3D plots of Pmiq)= " the probability of m paraboson charged-pairs 

-|- q positive parabosons" versus < n > and < >, the p = 1 curve is found to lie on the 

•i-H 

^ • relatively narrow 2D p = 2 surface. 
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1 Introduction: 

This paper is a theoretical study of multiphcity signatures in coherent production of charged-pairs 
of parabosons [1-3] of order p = 2 in comparison with those arising in the case of ordinary bosons, 
p — 1. The investigation is partially motivated by phenomenological coherence analyses, circa 
1970, of inelastic tt'^tt" pair production from fixed targets with laboratory kinetic-energies up to 
27 GeV by CP. Wang [4], and by Horn and Silver [5]. The present analysis is possible because 
the conserved-charge boson coherent states of [5, 6] are analogous to ones recently constructed for 
order p = 2 parabosons [7]. One physical consequence of "order p = 2" is that two or less such 
bosons can occupy a totally anti-symmetric state. Appendix 1 contains a brief review of p = 2 
paraboson statistics. 

Both the p = 2 and p = 1 models considered in this paper are notably "statistical" and crudely 
"unrealistic" in being free of kinematic and dynamical parameters associated with physically 
important production quantities such as the distributions of available energy/momentum, the 
kinematic size and other characteristics of the production region, resonances, masses and other 
conserved quantum numbers (e.g. isotopic spin in the p — 1 case [8]), because only the ^7(1)- 
charge conservation constraint has been imposed. Consequently, unlike what some readers would 
expect the situation to be from consideration of other particle-production models, in the figures 
the relative sizes of different aspects of the peaks, and of other structures, of the surfaces/curves 
cannot be adjusted without additional assumptions and complications. So despite the authors 
efforts, an interested viewer must sometimes force oneself to "look" in examination of some aspects 
of the figures. In the p = 1 case, more realistic models/arguments were made [8] to investigate 
such issues as independent-emission and coherence in 7r"^7r~ pair production but these analyses 
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are by now, of course, extremely primitive versus the current very sophisticated computational 
treatments [9, 10] of multi-particle production in high-energy collisions. For instance, even in 
electron-positron collisions a wide variety of coherence phenomena such as stage-one and stage- two 
spin-correlations in fermion-antifermion pair production, Bose-Einstein interference effects, and 
color-reconnection effects are used in various computational frameworks to describe and investigate 
particle production at LEP. Despite the lack of inclusion of other physics that would be required in 
a realistic application, we think that the focus in this paper on a comparison of these parameter-free 
statistical models is interesting and instructive with respect to p — 2 versus p — 1 parastatistics, 
and to the associated coherences involved in the production of conserved-charge paraboson pairs. 

Prom the then available data, in [4] Wang reported a regularity when, for various laboratory 
kinetic energies of the incoming primary particle, the relative frequency of events Pmio) was 
plotted against the mean multiplicty < n > for charged vr+vr^ pairs for pp, n^p, pn, n^n, and 
nn collisions. The data set was a compilation of about 50 inelastic production experiments with 
incoming primary kinetic energies extending from below production threshold to 27GeV. Wang 
tried to fit the regularity with two Poisson-type distributions. In [5], Horn and Silver constructed 
instead a simple parameter-free statistical model and found that, ignoring fluctuations in the data 
due to the specific structure of the various channels, the model did agree [5, 7] with the universal 
trends of the experimental regularity reported by Wang. One construction of the Horn-Silver 
model is/was from the p = 1 conserved-charge coherent states. Versus the inelastic reactions' 
multiplicity data, in [5] Horn and Silver argued for a statistical treatment of the gross features 
because (i) momentum conservation should be a weak constraint since the emitted pions occupy a 
small part of the available phase space, (ii) total isospin conservation on the distribution of charged 
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pions should also be weak since neutral pions are summed over, and so (iii) charge conservation 
remains as the important constraint. They derived the Pm{Q) distribution discussed below. Their 
Pm{0) distribution, i.e. the percentage of events with "m" 7r"'"7r~ pairs and no extra tt^'s , had 
been considered earlier by Kastrup [11]. 

At the present time, the only way known to construct the analogous model in the p — 2 case 
is to use conserved-charge coherent states. An alternative formulation with respect to either the 
statistical aspects or the coherence aspects is not currently available. So we now briefly digress to 
discuss what are these conserved-charge coherent states and then compare them with their p — 1 
counterpart: 

A charged-paraboson pair in order p — 2 consists of one A quanta of charge '-I-1' and one B 
quanta of charge '-1'. The Hermitian charge operator is defined by 

Q^N,-N, (1) 

where Na,b are the parabose number operators, see Appendix 1. Although Q does not commute 
with a or b, and although the paraboson pair operators ab ^ ba, since 

[Q, ab] = 0, [Q, ba] = 0, [ab, ba] = (2) 

the p = 2 coherent state can be defined as simultaneous eigenstate of Q, ab, and ba: 

Q\q, z, z' >— q\q, z, z' >, ab\q, z, z' >— z\q, z, z' >, ba\q, z, z' >— z'\q, z, z' > (3) 

Note that unlike in the p — 1 case (see below) , here two complex numbers z and z' arise because 
ab and ba are fundamentally distinct operators. Consequently, in the following multiplicity consid- 
erations, two non-negative parameters occur which are the moduli of these two complex numbers, 
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u = \z\ and v = \z'\. The explicit expressions for \q,z,z' > which involve the modified Bessel 
functions I[i]{u) and I^q+i^{v) are given in [7] for g > and q < 0; note [|] = "integer part of 
I" . For g > 0, the Mth-moment of the mean multiplicity for a p = 2 charged-paraboson pair is 

< >=< q, z, z'\{Ni))^\q, z, z' >. For q fixed, the percentage of events with m such pairs, Pm{q), 
is the square of the moduli of the expansion coefficients of \q,z,z' > in terms of the two-mode 
parabose number Fock states ( see [7] ). The explicit multiplicity formulas corresponding to the 
figures shown in this paper can be found below in Section 2. 

The simpler p — 1 conserved-charge coherent state \q,^ > satisfies ^ >= q\q,^ > and 
ct^l?)^ >= Cl?)^ >• Since a6 = 6a for p = 1, \q,^ > involves a single complex number ^, and so 
the Horn-Silver multiplicity analysis involved a single, non-negative real parameter x = |^|. 

Because the p — 2 construction involves two non- negative real parameters u and v, we consider 
a 3D plot of Pm{q) = Z{u, v) versus the quarter plane defined by < n >=< Ub >= X{u, v) and 

< >=< {ribY >— Y{u, v). This is the obvious generahzation of Wang's plot and it is the focus 
of Section 2 and its figures. Section 3 contains some additional multiplicity formulas which are 
similar in structure for p — 1,2. Section 4 contains a brief discussion of some of what has/hasn't 
been learned from this analysis of coherent production of charged-pairs of parabosons of order 
p = 2. Appendix 2 and its figures are focused on how different shaped {u, v} charts, or coordinate 
patches, appear when mapped onto the two-dimensional surface. 



2 Pmisi) versus < n >, < >: 

The analysis in this section is straight-forward because there are simple explicit mathematical 
formulas involving modified Bessel functions and because "Mathematica" [12] notebooks are con- 



venient for displaying the parametric 1-parameter p — 1 "curve" and 2-parameter p — 2 "varying- 
width, folded ribbon." We find that a concise "curve/ribbon" terminology is useful and appropriate 
for describing, for example, Figs. la,b, which are for m = 2, and q — I. On a 2D side of the 3D 
display, the p = 2 ribbon consequently appears projected as a relatively narrow, two-dimensional 
region which will be denoted as a "band" ; see Fig. lb which shows the front XZ side of Fig. la. 
In spite of their being from the more complex u v asymmetric case, these two lead figures have 
been chosen in order to ab initio discuss some details and caveats. 

If the reader initially finds Figs. la,b and their description too complex, the reader might go 
directly to Figs. 2 and 3a,b which are for the simpler u ^ v symmetric case. 

In Figs. la,b, and in the remaining figures in this paper, the solid line is the p — 1 curve. In 
these two figures, this p — 1 curve is also the {0, v} line on the p — 2 ribbon. Near the peak, 
the "fold" is the bottom edge of the ribbon. As shown, the p = 2 ribbon consists of open-circles 
for the non-folded u > v region, and of solid-circles for the folded u < v region. The upper 
edge of the ribbon is the {u, 0} set of points. Slightly to the right of the peak, one can see from 
the solid-circles that each hne of dots traveling leftward down the page, bends under (or "over" , 
whichever as the viewer prefers) the fold to reach the {0, v} line on the ribbon. 

For q odd, we assume it is a "complete fold". The words "complete fold" mean that one part 
of the {u,v} ribbon lies exactly on top of another part, and so there is no difference in going 
"under" or "over" such a fold. 

Due to the very limited analytic results, mental/visual care is generally advisable: (i) Only for 
the cases of g = and for q — 1 is it proven that the p — 1 curve always lies on the p — 2 two 
dimensional surface. In these two cases the proof appears un-instructive for then the parametric- 
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solution u — 0,v — 2x { and for g = also u — 2x, v — ) avoids confronting different ratios of 
modified Bessel functions. For an arbitrary q value, from the asymptotic limits for the modified 
Bessel functions, it follows that as x — > 0, the curve and as v — > 0, the ribbon must begin at 
the same point: the origin for m 7^ and Po{q) = 1 for m = 0. Similarly, for large parameters for 
arbitrary q, the curve and ribbon will also be at the same 3D point with Pm{<l) — > 0, if x = \{u-\-v). 
(ii) Only for q even is it known that there is a complete fold of the ribbon because then Pm{(l), 
< n >, and < v? > are each u ^ v symmetric. For q even, the fold is the line u = v, c.f. the 
simpler Figs. 2 and 3a,b. 

The "line of dots traveling leftward down the page" discussed above for the peak region in 
Figs. la,b are a set of {u, v} values from a unit-negative-slope diagonal in the {u, v} domain, see 
last figure of this paper and its discussion Fig. 8a,b. 

Second, the reader should be aware that the "dots" displayed to show the two-dimensional 
surface correspond to specific {u,v} parameter values and that the associated ui,U2, vi,V2, ... 
values are equally spaced, c.f. Fig. 5a below where some {u, v} coordinate values are shown. In 
some figures, U2 — Ui is not equal to V2 — Vi. We find that a careful usage of a few un-connected 
"dots" does not mislead in displaying these two-dimensional surfaces. The "dots" should not be 
confused with a random generation of data points, such as in a scatter plot. Note that the u and 
V parameters do not map into an orthogonal coordinate chart on the ribbon. In general, at each 
point on the ribbon the curvature is non-zero. 

Third, in consideration of Figs. la,b, as well as the other figures in this paper, it should be 
noted that, even near the peak, the width of the ribbon is narrow versus, for instance, the half- 
width of the p = 1 peak. This fact, besides the folding and the role of the u and v parameters. 
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would be an important issue in attempting empirically to distinguish coherent production pairs 
oi p — 2 parabosons versus those of ordinary p = 1 bosons. On the other hand, the ribbon width 
is indeed non-zero and to the eye there is generally a systematic above, or below, displacement of 
the visual center of the p = 2 "band" from the p = 1 curve as in Figs. la,b. 

For the m — 2, and q — 1 case shown in Figs. la,b the formulas for the p—1 curve are 

and for the p — 2 ribbon are 

2 Iq[u) Ii{v) 

< >(^) = < n >(^) +]{^ + + (7) 

4' Io(u) Io{u)h(v) h(v) 

(2) _ 3u^v + Gu^^ + 

384/o(^)/,(.) 

As in (4-8), a superscript "(1) or (2)", for p = 1 or p = 2, can be respectively written on < n >, 
< >, and Pmio) when needed to avoid confusion. Usually one knows the p order from the 
working context, and so these superscripts can often be suppressed. 

Before discussing other specific m and q cases, we consider the parametric formulas for an 
arbitrary curve and ribbon. 

2.1 p = \ curve: 

The p—1 curve {X{x),Y[x), is parametrized by the real non-negative parameter x . For q 

> 0, in terms of modified Bessel functions, Iq{2x), 

^ <">=^ 
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Y{x) = <v? >=< n > -F 



Iqi^x) 



(10) 



Equivalently, < >— x^{l — ^-^fj^^) follows due to the recursion relation for J^'s. The proba- 





bility for " m boson pairs + q positive bosons " is 



X' 



.2m+q 



Z{x) 



Pm{q) 



Iq{2x)m\{m + q)\ 



(11) 



For example, in the analysis of "ideal" data for the production of purely multi-pion final states, 
Pmio) would be the probability for the production of "m + g" tt+'s and "g" 7r~'s. For p = 1, it is 
instructive to consider the statistically-fundamental 



that it is via normalization that the modified Bessel function appears in development of the 
simple idea of charge-conservation-constrained Poisson distributions for independent tt^ and vr" 
production [5, 8]. 

For multi-pion final states Wang [1], and later Horn and Silver [2] in consideration of (12) as a 
model, plotted Pmil) versus < n > for small values of m and q. This is the XZ plane, or "front" 
projection of the 3D figures considered in this paper. 

2.2 p = 2 surface: 

The p— 2 two-dimensional surface {X{u^ v),Y{u, v), Z{u^ v)) is parametrized by the two real non- 
negative parameters u and v. As for the p = 1 curve discussed above, only the Z-coordinate is 
m dependent. 




,2m 



(12) 



with an m-independent normalization constant (Ng) 



-2 



X ^/g(2x) because in this way one sees 
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{q Even): The case of q even and g > is symmetric in u v with Z{u,v) — Z{v,u), see 
below, and the u v symmetric 

I f Uls+2{u) Vls+2{v)\ 

YM . < n >(-) +1 ( + 2. /f/^^ ^f/^^ + (14) 

{q Odd): The case of q odd and > is asymmetric in tt <-> v with Z{u,v) ^ Z{v,u), and 
the u ^ V asymmetric 

1 / Ul3+l{u) Vl3+3{v)\ 

X{u,v) ^ <^>^'^=^ T^ + T^ ^1^) 

1 (u^Is±3{u) Is±i{u)I^{v) V^I,_±5{v)\ 

Y{u,v) = < >(^)^< n >(^) +- ^ + r / ^ r / ^ + r %a ^6 

4 1 /a^(u) (u) /a+i (w) /2+i(t') / 



2 2 2 



For g > 0, the probabihty for " m paraboson pairs + q positive parabosons " is 

m+l 

Ziu,v) ^ Pi^H?) = (A^f f E (17) 

1=1 

(jV«)-2 = (18) 

P . = (19) 

where in Pq^rn;i 

r . r-'-)'""^-'^ !-(-)'] (20) 

. , + (-)'] pi) 

In paraboson statistics, it is sometimes convenient to use [x] — "integer part of x", and [x\\ — 
[x\[x — 1] • ■ ■ 1, [0]! = 1. Note that this [x] symbol occurs in (18-21). The use of the symbol [x] 
enables a compactification of several analytic expressions and its use is convenient in computational 
notebooks. However, it is sometimes instructive, as in (13-16), to write out the various cases 

9 



by letting q — 2(5, or 2(5 + 1; Q = 0, 1, • • •; or similarly for m or i. Throughout this paper, 
the appearance of "square braces" always denotes such an integer truncation. Note that the 
summation over "i" in (17) is also due to the fact that ah ^ ha [7]. The compact expressions are 



2 



Since from (17) the p — 2 case non-trivially involves two parameters, u and v, the figures in 
this paper are 3D ones with P^\q) = plotted versus the two-dimensional < n >(^^ = X^^ 
and < = Y^^^ plane. Only the first octant is used. 

2.3 Other m and q cases: 

For the case m — 0, q — the parametric expressions are very simple: 

<n^>-=<n>-.^ (-) 
^°'(°)-io(W ^^'^ 

< >(^)=< n >(^) +i(^ + + !^) (25) 

Io{u)Io{v) 

This is the simplest u ^ v symmetric case and the "right" FZ-side is shown in Fig. 2. As 
previously mentioned, in the symmetric case, the fold is along u = v which is shown by the dashed 
fine in Fig. 2. The two edges of the original ribbon, {0, v} and {u, 0} are respectively mapped by 
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V — 2x, u — 2x, into the p — 1 curve. Only open-circle points from the u > v region are shown. 
All m — curves/ribbons drop down from Po{q) — 1. 

Figs. 3a,b, respectively show the 3D plot and the right YZ-side for another symmetric case 
m = 1, q = 0. Note that, unlike in the preceding Fig. 2, in this case the p = 1 curve now lies on 
the top of the ribbon until there is an almost 180° twist, after the peak where the p — 1 curve is 
on bottom. In Fig. 3b the FZ-side quite clearly shows this twist. The parametric expressions are 
respectively A«(0) = ^ , Pi'\o) = 

In appendix 2, there is a discussion of the {u, v} coordinate charts associated with the asym- 
metric case m = 1, g = 1 shown in Figs. 4a - 4d. The parametric equations are respectively 

Pi '{1) = 2i^(2x) ■> -^1 (-*-) ^ 8/o("Ki(f) • [ '^^^ omitted figures for the asymmetric case m = 0, g = 1 
similarly show the fold occurring in the u < v region. ] Fig. 4a is the 3D display and Fig. 4b 
is the projection onto the front XZ plane. Fig. 4c is from a different 3D viewpoint. It provides 
a view back into the origin, with a YX floor below and a YX ceiling above, in a "left-to- right, 
front axes labeling" . This figure shows how nearly vertical the ribbon is. Lastly, Fig. 4d shows a 
close-up of the m = 1, g = 1 peak. This figure illustrates the fold. In it the upper sohd line is the 
p — 1 curve and the lower solid line is the u — v curve which separates the open-circles {u > v) 
from the solid-circles {u < v). 

The different shaped {u, v} charts of Appendix 2 can be used in investigating the folding which 
occurs in a transverse crossing of the ribbon. Similarly, the lines u = 0, also v = for g-odd, 
u — V, and for g > 2 also the independent p — 1 curve, can be used. Both techniques are useful in 
studying what occurs in going out along the ribbon. In proceeding from the origin, or for m = 
from Po{q) = 1, there are twist (s) and a peak (absent for m = 0). If as in Fig. 4c, the ribbon is 
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viewed back in from the high u and v parametric coordinates, then in moving out from the origin 
but viewing back towards the origin, the almost 180° twisting of the u > v edge from the peak 
onward is counter-clockwise for (m, q) — (1, 0), (1, 1), (2, 1). For (2, 1), (2, 2), (2, 3) there is also a 
twist near the origin. For the peak in (1, 1), c.f. Fig. 4a, the Z coordinate values are: for the 
V = top edge, Z(2.65,0) = 0.475; for the p = 1 curve, Z(0,3.25) = 0.436; for the u ^ v curve 
Z(1.85, 1.85) = 0.42; and for the fold Z(1.5, 2.25) = 0.41. 



3 Additional Multiplicity Formulas: 

In the p = 2 case, the M-th moment of the mean multiplicity 

< < {N,nD;,^r{N,)-' (27) 

where IDi, — + v-^)^ Da — Df,-\- q. It follows that 

<n>^<nb>^ {NgfSi{u,v), (28) 

< >=< nl >= iN,)\S2 + Si), (29) 

< nl >= {N.fiSs + 3S2 + S1), (30) 

< nt >= {Ng)\S4 + 6S3 + 7S2 + Si), (31) 
< nl {Ng)^{S5 + lOS^ + 25S3 + I5S2 + Si), ... (32) 

where, with a binomial coefficient, 

Sm{u,v) = cE(^)(|)^-'(|)'/[.±H(|£^](«)/[^](^^) (33) 

c ^ 2[f]+[^l«-[fli;-[^l (34) 
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The function Sm+i{u,v) is generated by (D^ — M)Sm{u, v) — Sm+i{u, v). The expression for 
{N^y is (18) above. 

Note that as in the p—1 case, 

<-r>-E('^)^*<-r*> (35) 

When q = 0, < n*^ >=< >, and for fixed q, < Ua >=< rib > +<!■ 

For (q— — \q\) < 0, there are the following relations to the g > expressions: 

< nf >_|,|=< nf >|,| (36) 

in (27-32) and exchange u ^ v, 

< >_|,|=< nf >,„ (37) 

in (27, 35) and exchange u v; and Pm(|g|) = Pm+\g\{—\q\) . 
The analogous formulas for the p—1 case are 

< < >^^^= im\^ir{N(^Y' (38) 

where D^''^ = \{x£),D^a^ = D^^^ + q. So < n >(^)=< rib >^^^= {N^^^fsi{x) , 
< nl iN^^^V{s2 + si), ... where 

sm{x) = x^-^I.+Mix) (39) 

The function sm+i{x) is generated by (-D^^'* — M)sm{x) = sm+i{x). The expression for (A^^^^)^ is 
given after (12) above. 

For {q = -\q\) < for p = 1, < >-\q\^< >|,| , < >-\q\^< >|,|, and 

PM)=Pm+MHQ\)- 
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4 Discussion: 

The 2D surface occurs in the 3D plots of the relative probabihties Pm{q) versus < n >, < > 
as required by the coherences embodied in the equations in (3) for the p = 2 conserved- charge 
coherent states. Unlike the analogous p — 1 case, where such a model was motivated from 
a regularity in tt+tt" data, it is necessary in the case of the production of the charged-pairs of 
parabosons of order p = 2 to assume an analogous coherent production mechanism and to assume a 
reasonable, but definite, treatment of the ab and ba operators. Neither of these two assumptions 
might be true if paraboson pairs are found to be produced in nature in high energy physics 
collisions, e.g. in central-diffactive-exchange experiments [14], or in some area of contemporary 
astrophysics/cosmology, e.g. in dark matter detectors [15] or by very high energy cosmic rays 
[16]. In the infra-red domains of QED and of QCD, coherent-state/degenerate-state coherence 
is a well-known phenomena. On the other hand in QCD, partonic jets dominate the very high 
energy hadronic production processes. In the case of the lower energy tt+tt" production, as was 
emphasized by Horn- Silver, and others, the situation is complex and there is much more physics 
and phenomenological structure than that incorporated in the reference p — 1 model considered in 
this paper. In the approximation in which one neglects such additional physics, this p = 1 model 
is parameter free because the x parameter is effectively replaced by < n >, the mean number 
of final charged pairs. Similarly, the analogous p — 2 model is parameter free for the u and v 
parameters are effectively replaced by < n > and < ri^ >. At the present time, the existence and 
relative importance of other conserved quantum numbers, of resonances or other phenomenological 
interaction effects, etc. is completely unknown for a coherent production of pairs of parabosons 
of order 2. Nevertheless, with respect to experiments [14-16], one conclusion from this paper is 
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that from a parameter-free statistical model, one would expect a signature of relatively narrow 
bands, due to projection of varying-width folded ribbons, will be present in 2D analysis of pair 
multiplicities from coherent production of parabosons of order 2. 

In the case of p = 2 parabosons, there will be a second kinematic variable [4] and multiplicity 
data will not scale in terms of a single- variable curve [11]. For instance, if the variable for the 
p — 1 curve is X ~ ^{u + v) , then for p = 2 parabosons there will also be a dependence on 
y ^{u — v). In Fig. 8a, with x ~ ^{u + v) fixed, the sensitivity to the second variable 
y ~ |(m — is very striking and significantly greater than the naive width of the ribbon. In the 
region of the peak. Figs, lb and 4b show that this signature is considerably enhanced when there 
is an extra positive paraboson, versus the case of only charged-paraboson-pairs. Fig. 3b. While 
an explicit dynamical production model is required to theoretically investigate the sensitivity with 
respect to specific kinematic variables, in independent-particle-emission models the total energy 
in the emitted particles is monotonically related to the intensity strength of the source. This 
suggests that two simple and useful kinematic variables are the sum, Etotai = total + E~ total, 
and difference E'^fotai ~ E'totah of the total emitted Q — 1 and Q — —1 energies. 

The "identical particle" defining parabose tri-linear relations leave the two pair operators, ab 
and ba, fundamentally distinct which (3) and present analysis maintain by the use of the two 
distinct u and v parameters. In the p — 1 case, the parameter can be interpreted as 

a product of the intensity strengths of the tt"*" and 7r~ sources. In the p = 2 case, u and v can 
be interpreted as the intensity strengths of the "a6" and "6a" sources. Their difference is due to 
the two distinct orderings of the a and b operators. The physical observables are not always u <-> 
V symmetric. In particular, the presence of a g odd total charge produces an asymmetry in the 
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folding of the 2D surface; for instance, see Eqs. (15,16) and (8). 

There is overall <-> B~ symmetry {U{1) charge symmetry ): The "Probability for (m+ 
74+'s and (m) S~'s " equals "Probability for (m) A+'s and (m + |g|) 5~'s ", because Pj„(|g|) = 
Pm+\q\{-\q\) and < >_|g|=< n^^ >\q\ for M =integer. 

The U{1) charge might be a "hidden conserved-charge" such as to yield pure pair production, 
i.e. only final state events with g = would occur such as in the production of strange particles 
via the strong interactions. In this case, or by a designed selection of only pair final state events, 
the p — 2 versus p — 1 results of this paper can still be used with Z^^=o-fm(0) = 1. This is 
the simpler symmetric case and because of g = there would be more analytic control in such an 
analysis. 
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Appendix 1: Parabosons of order 2 

In local, relativistic quantum field theory, identical particles obey either (i) parabose and 
parafermi statistics [1-3] for which the number of particles in an antisymmetric or a symmetric 
state, respectively, cannot exceed a given integer p, or (ii) for two space dimensions, infinite 
statistics based on the braid group [13]. In parabose statistics, instead of bilinear, there are 
fundamental trilinear commutation relations: 

[ofc, {aj, a„}] = 26kiam, [ak, {a}, a}^}] = 25kia}^ + 24^a]', 

(40) 

[ak,{ai,am}]^0, (A;,/,m = l,2) 
where [C, D] = CD — DC, and {C, D} = CD + DC. In the case of only two kinds of parabosons, 

there are some simple commutation relations between "A" paraboson operators and the "B" 

paraboson operators: letting a = ai, b = a2, 

[a, = 0, [b, a^] = 0, [a^ 6'] = 0, [6^ a"] - (41) 

plus the hermitian conjugate relations. Order p — 2 is simpler because there is the "self-contained 
set" of 3 relations [2] 

ttmaial - o-k^i^^m = 2 Ski a-m. (42) 

akA^^m - ama\ak = 2 4; - 2 5im ak (43) 

akdiam - dmaio-k = (44) 
The parabose number operators for p = 2 order are defined by 

Na = \{a\ a} - 1, TV, = \{b\ b} - 1 (45) 

In [7], the state- vector space for two-mode parabosons and the order p — 2 conserved-charged 
parabose coherent states were constructed. 
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Production and decay selection rules [2, 3] exclude known particles from obeying other [1] than 
the usual boson and fermion statistics, i.e. order p = 1 parastatistics. 
Appendix 2: {w, v} Chcirts on 2D surface 

In a determination of the statistical uncertainties for an application of the parametric formulas 
in Section 2 in comparison with actual experimental data, the structure(s) of the folded, partially- 
twisted ribbons versus p — 1 curves will enter and have to be correctly treated. For this reason, 
and for possible more abstract theoretical/mathematical use by other readers, the following details 
about how different shaped {u, v} coordinate patches appear when mapped onto the 2D surface 
should be of interest. The lines u = 0, also v = for g-odd, and u = v provide useful standard 
references in examining what is taking place on each ribbon. In addition, for g > 2 the p — 1 
curve is a different reference since it no longer is the u — line. 

The m = 1, q = 1 case which is u ^ v asymmetric is used in this appendix. This case is both 
generic and visually simple. Figs. 4a-4d were earher displays of this case. The discussion begins 
with the non-folded u > v region: 

Figs. 5a,b show how the u > v "strip" from u = 2 to u = 3 is mapped onto the peak region of 
the ribbon. In 3D Fig. 5a, the solid line that lies on the bottom of the "patch" on the ribbon is 
the u = V set of points. The open-circle points on the top of the patch are the = set of points. 
The points on the origin side in Fig. 5a are for the vertical u — 2 line in Fig. 5b. The points on 
the patch farthest from the origin, which appear to have a nearly constant negative slope, are for 
the other vertical line in Fig. 5b, i.e. for u — 3. Here there have been 5 sets of vertical points 
mapped from Fig. 5b onto the 5a patch on the ribbon. 

Figs. 6a,b similarly show the differences when the u > v points arise from three more widely 
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separated constant-w vertical lines with respectively u ~ 2,3,4. To the left-side of the peak. 
It ~ 2, the spacing between these constant-w-value points is rather small, whereas on the right of 
the peak, where it ~ 4, the spacing between such points is significantly greater. The points are 
equally spaced in v so there are more points taken at w ~ 4 than at m ~ 2. 

Horizontal "strips" can similarly be used to study the folded u <v region. However, a "ring" 
region or a "unit-negative-slope strip" region can be used which includes both u > v and u < v 
domains. Figs. 7a,b show what occurs when points in a "ring" of radii r< = 3 and r> = 4 are 
mapped onto the ribbon. In 7a, the open-circle-dots are again from the simpler u > v region, 
whereas the solid-dots are from the u < v region where the fold occurs. Fig. 7c shows a close-up 
where the u — points {p — 1 curve) are on the upper solid curve; the u — v points which divide 
the ring in half are on the lower solid curve. [ If instead, the points are from pie-slice-like sectors, 
the associated lines on the ribbon then extend indefinitely on out the ribbon. ] 

Lastly, Fig. 8a,b show what happens when points are taken from a somewhat similar, but 
sometimes more useful, unit-negative-slope strip in the u, v quadrant. In Fig. 8a the solid line is 
the u — v curve. The "fine of dots travehng leftward down the page" discussed in the text for the 
peak region in Figs. la,b are a set of {u, v} values from such a unit-negative-slope diagonal. If 
instead, the points are taken from the orthogonal unit-positive-slope strips, the associated lines 
on the ribbon will start on a {0, v} or {u, 0} line, depending on which side oi u — v the strip is 
from, and then the line will extend indefinitely on out the ribbon. 



19 



References 

[1] H.S. Green, Phys.Rev. 90, 270 (1953); D.V. Volkov, Sov. Phys.-JETP, 11, 375 (1960); and 
c.f. W. Pauli, Handbuch der Physik Vol. 24, p. 113 of Sec. 14 (1934); p. 128-9 in A Scientific 
Biography of W. Pauli, CP. Enz (2002) (Oxford U. Press). 

[2] O.W. Greenberg and A.M.L. Messiah, Phys.Rev. 136, B248 (1964); ihid., 138, B1155 (1965). 

[3] Y. Ohnuki and S. Kamefuchi, Quantum Field Theory and Parastatistics (1982) (Univ. Press 
of Tokyo). 

[4] CP. Wang, Phys.Rev. 180, 1463 (1969). Wang plotted P^{q) for g = versus the equivalent, 
mean number of charged prongs < Hs >— 2 < n > +2 = < Uc > (Horn-Silver). Figure la in 
Wang's paper shows the laboratory kinetic energies of the data points. 

[5] D. Horn and R. Silver, Phys.Rev. D2, 2082 (1970). 

[6] D. Bhaumik, K. Bhaumik, and B. Dutta-Roy, J. Phys. A: 9, 1507 (1976); J.R. Klauder and 
B.-S. Skagerstam, p. 43-48, Coherent States{19S5) (World Sci. Pub. Co., Singapore); and 
see G.S. Agarwal, J. Opt. Soc. Am. B5, 1940 (1988). 

[7] S. Jing and CA. Nelson, J. Phys. A: 32, 4131 (1999). 

[8] See, and references therein, D. Horn and R. Silver, Annals of Phys. (NY) 66, 509 (1971); D. 
Horn, Physics Rep. 4, 1 (1972); and D. Horn and F. Zachariasen, Hadron Physics at Very 
High Energies (1973) (W.A. Benjamin, Reading, MA). 

[9] T. Sjostrand and M. Bengtssson, Comp. Phys. Comm. 43, 367 (1987); G. Marchesini and 
B.R. Webber, ihid. 67, 465 (1992); S. Jadach, Z. Was, R. Decker, and J.H. Kuhn, ihid. 76, 

20 



361 (1993); S. Jadach, B.F.L. Ward, and Z. Was, ibid. 124, 233 (2000); B.F.L. Ward and 
S. Jadach, Acta Phys. Polon. B33, 1543 (2002); L. Lonnblad, Comp. Phys. Comm. 71, 15 
(1992); and T. Sjostrand, P. Eden, C. Priberg, L. Lonnblad, G. Miu, S. Mrenna and E. 
Norrbin, ihid. 135, 238 (2001). 

[10] See, and references therein, Perturbative Quantum Chromodynamics, ed. A.H. Mueller (1989) 
(World Sci. Pub. Co., Singapore); Yu L. Dokshitzer, V.A. Khoze, A.H. Mueller and S.I. 
Troyan, Basics of Perturbative QCD, ed. J. Tran Thanh Van (1991) (Gif-sur-Yvette, Editions 
Frontieres); J.C. Collins and D.E. Soper, Ann. Rev. Nucl. Part. Sci. 37, 383 (1987); R.K. 
EUis, W.J. Stirhng and B.R. Webber, QCD and Collider Physics (1996) (Cambridge Uni. 
Press); B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rep. 97, 31 
(1983); B. Andersson, The Lund Model (1997) (Cambridge U. Press) ; L.L. Frankfurt and 
M.I. Strikman, Phys. Rep. 160, 235 (1988); L.N. Lipatov, ibid. 286, 131 (1997); I. Dremin 
and J.W. Gary, ibid. 349, 301 (2001); S. Catani, F. Krauss, R. Kuhn, and B.R. Webber, 
JHEP 0111, 63 (2001); C.F. Berger, T. Kucs and G. Sterman, Phys. Rev. D68, 014012 
(2003); A. Hebecker, Phys. Rep. 331, 1 (2000); and E.A. De Wolf, J. Phys. G: 28, 1023 
(2002). 

[11] H.A. Kastrup, Nuc. Phys. Bl, 309 (1967); H. Gemmel and H.A. Kastrup, Z. Phys. 229, 
321(1969); Nuc. Phys. B14, 566(1969). 

[12] Mathematica, Version 4.2, Wolfram Research, Inc. (Champaign, IL). 



21 



[13] Fractional Statistics and Anyon Superconductivity, ed. F. Wilczek (1990) (World Sci. Pub. 
Co., Singapore); and A. Khare, Fractional Statistics and Quantum Theory (1997) (World 
Sci. Pub. Co., Singapore). 

[14] M.G. Albrow and A. Rostovtsev, hep-pli/0009336; V. Khoze, A.D. Martin, and M. Ryskin, 
Eur. Phys. J. C14, 525(2000); hep-ph/0207313; A. De Roeck, V. Khoze, A.D. Martin, et.al, 
Eur. Phys. J. C25, 391(2002); CDF CoUab., D. Acosta, et.al., hep-ex/0311032; TOTEM 
CoUab., http://totem.web.cern.ch/Totem. 

[15] DAMA CoUab., R. Bernabei, et.al, Riv. N. Cim. 26, 1 (2003); CDMS CoUab., R. Abusaidi, 
et.al., Phys. Rev. Lett. 84, 5699(2000); EDELWEISS CoUab., A. Benoit, et.al., Phys. Lett. 
B513, 15(2001); IGEX CoUab., H. Morales, et.al., ibid., B489, 268(2000); HDMS CoUab., 
H.V. Klapdor-Kleingrothaus, et.al, Nuc. Phys. Proc. Supp. 110, 364(2002); ZEPLIN Col- 
lab., J. Dawson, et.al, 110, 109(2002); SIMPLE CoUab., J.I. Collar, et.al. New J. Phys. 2, 
14(2000); PICASSO CoUab., N. Boukhira, et.al, Nuc. Phys. Proc. Supp. 110, 103(2002). 

[16] AGASA CoUab., M. Takeda, et.al, Astropar. Phys. 19, 447 (2003); Hi Res Col- 
lab., D.R. Berman, et.al, Nuc. Phys. Proc. Supp. 117, 106 (2003); AUGER CoUab., 
http://www.auger.org/; and c.l V.A. Karmanov and A.E. Kurdyavtsev, hep-ph/0207321. 

Figure Captions 

Fig. la: The probabiUty -P2(l) of "2 charge-paraboson-pairs plus 1 positive paraboson" versus 
< n >, < >. This is the more complex u ^ v asymmetric case. The solid line is the p = 1 
curve. Near the peak, the fold is the bottom edge of the ribbon. The open-circles are for the 
non-folded u> v region, and the solid-circles for the folded u < v region, see discussion in Sec. 2. 
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Fig. lb: The front projection, or XZ-side, of Fig. la, so ^(l) versus < tt. >. 

Fig. 2: The right-side projection, or FZ-side, for the probabihty -Po(O) of "zero charge- 
paraboson-pairs plus zero additional single charged parabosons" versus < ri^ >. Since q is even, 
this is the simpler u ^ v symmetric case. Only open-circles from the u > v region are shown. 
The fold is the dashed hne u — v. 

Fig. 3a: The probability -Pi(O) of "1 charge-paraboson-pair plus zero additional charged para- 
bosons versus < n >, < >. This is also the simpler u ^ v symmetric case so only open-circles 
from the u>v region need to be shown. The sohd line is the p—1 curve. The fold is the dashed 
line u = V. 

Fig. 3b: The right projection, or FZ-side, of Fig. 3a, so Pi(0) versus < >. 
Fig. 4a: The probability -Pi(l) of "1 charge-paraboson-pairs plus 1 positive paraboson" versus 
< n >, < >. This case is u ^ v asymmetric. 

Fig. 4b: The front projection, or XZ-side, of Fig. 4a, so -Pl(I) versus < n >. 
Fig. 4c: A 3D view "back into the origin" of Fig. 4a. 

Fig. 4d: A close-up view of the peak of Fig. 4a. Appendix 2 and its figures discuss various 
{u, v} coordinate charts for this asymmetric case. 

Fig. 5a: For Pi(l) of Figs. 4, the image of the u > v "strip" from w = 2 to « = 3 shown in 
Fig. 5b. 

Fig. 5b: The domain u> v "strip" from u = 2 to u = 3 for Fig. 5a. 
Fig. 6a: For Pi(l), the image of the three vertical hues in Fig. 6b. 

Fig. 6b: The three vertical hues u ~ 2, 3, 4 which are mapped to the 2D surface in Fig. 6a. 
Fig. 7a: For -Pi(l), the image of the "ring" of points of Fig. 7b. 
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Fig. 7b: The domain of points in the "ring" of radii r< = 3 and r> = 4 which are mapped 
into Fig. 7a. 

Fig. 7c: A close-up of Fig. 7a. 

Fig. 8a: For the image of the points in the unit-negative-slope domain in Fig. 8b. The 

{u, v} values shown are a subset from the unit-negative-slope line of points nearest the origin in 
Fig. 8b. Note the location of the end points {4, 0} and {0,4}. 

Fig. 8b: The unit-negative-slope domain whose points are mapped into Fig. 8a: 
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